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ABSTRACT 

We present a series of observations of the return of activity in main-belt comet 
238P/Read. Using data obtained in July and August 2010 when 238P appeared to be 
largely inactive, we find best-fit IAU phase function parameters of H = 19.05±0.05 mag, 
corresponding to a nucleus radius of r n ~ 0.4 km (assuming an albedo of pr = 0.05), 
and G = —0.03 ± 0.05. Observations from September 2010 onward show a clear rise in 
activity, causing both a notable change in visible morphology and increasing photometric 
excesses beyond what would be expected based on bare nucleus observations. By the end 
of the observing period reported on here, the dust mass in the coma shows indications 
of reaching a level comparable to that observed in 2005, but further observations are 
highly encouraged once 238P again becomes observable from Earth in late 2011 to 
confirm whether this level of activity is achieved, or if a notable decrease in activity 
strength compared to 2005 can be detected. Comet 238P is now the second main-belt 
comet (after 133P/Elst-Pizarro) observed to exhibit recurrent activity, providing strong 
corroboration for the conclusion that it is a true comet whose active episodes are driven 
by sublimation of volatile ice. 

Subject headings: comets: general — comets: individual (238P/Read) — minor planets, 
asteroids 



'Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California, and the National Aeronautics 
and Space Administration. The Observatory was made possible by the generous financial support of the W. M. 
Keck Foundation. Some data presented herein were also obtained at ESO facilities at La Silla under program ID 
085.C-0363(A). 
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1. INTRODUCTION 

Main-belt comets (MBCs) are a recently-identified population of solar system objects that 
exhibit cometary activity but are dynamically indistinguishable from main-belt asteroids (Hsieh & 
Jewitt 2006). Discovered on UT 2005 October 24 (Read et al. 2005), comet 238P/Read (hereafter, 
238P) was the second main-belt comet to be discovered. A physical study of the comet shortly 
following its discovery found that its activity was consistent with dust emission driven by the 
sublimation of recently excavated subsurface ice, and inconsistent with impact-driven dust emission, 
meaning that 238P is a true comet (Hsieh et al. 2009b), and not a disrupted asteroid (cf. Jewitt 
et al. 2010). That study further found a mass loss rate roughly an order of magnitude larger than 
that calculated for the first known MBC, 133P/Elst-Pizarro (Hsieh et al. 2004). This higher rate 
of dust emission may indicate that 238P's activity was triggered much more recently than that of 
133P, and therefore is driven by the sublimation of a much less depleted volatile supply (Hsieh et 
al. 2009b). 

Haghighipour (2009) found 238P to be dynamically unstable on a time scale of ~20 Myr. Its 
proximity to the Themis collisional family, however, suggests that it may have once been a stable 
Themis family member, but was destabilized (via a gradual increase of its orbital eccentricity) by 
the nearby 1:2 mean-motion resonance with Jupiter. As such, despite its current instability, 238P 
may be physically related, and therefore compositionally similar, to actual Themis family members 
and fellow MBCs 133P and 176P/LINEAR. 

2. OBSERVATIONS 

i?-band observations of 238P were made in photometric conditions between July and December 
2010 using the 4.2 m Southern Astrophysical Research (SOAR) telescope at Cerro Pachon and the 
3.54 m New Technology Telescope (NTT) at La Silla, both in Chile, and the 10 m Keck I and 
2.2 m University of Hawaii (UH) telescopes on Mauna Kea in Hawaii. Details of these observations 
are listed in Table 1. Observations with SOAR were made using the SOAR Optical Imager (SOI; 
Schwarz et al. 2004) which utilizes a mosaic of two E2V 2048x4096 pixel CCDs (0 / .'l54 pixel -1 using 
2x2 binning) behind standard Kron-Cousins broadband filters. NTT observations were made using 
the ESO Faint Object Spectrograph and Camera (EFOSC2; Buzzoni et al. 1984), which employs a 
2048 x 2048 pixel Loral/Lesser CCD (Of! 24 pixel -1 using 2x2 binning), behind Bessel broadband 
filters. Observations with Keck were made using the Low Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995) in imaging mode. LRIS employs a Tektronix 2048 x 2048 CCD with an image scale 
of / .'210 pixel -1 and Kron-Cousins filters. UH 2.2 m observations were made using a Tektronix 
2048 x 2048 pixel CCD (tf'219 pixel -1 ) behind Kron-Cousins filters. 

Standard bias subtraction and flat-field reduction were performed on all data. For SOAR, 
NTT, and UH 2.2 m data, flat fields were constructed from dithered images of the twilight sky, 
while images of the illuminated interior of the dome were used to construct flat field images for 
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Keck data. Photometry of Landolt (1992) standard stars and field stars was obtained by measuring 
net fluxes within circular apertures, with background sampled from surrounding circular annuli. 
Comet photometry was performed using circular apertures with radii selected to match the seeing 
conditions of each night, but to avoid the contaminating effects of the coma, background sky 
statistics were measured manually in regions of blank sky near, but not adjacent, to the object. 
Field stars in each comet image were also measured to assess and correct for extinction variation 
during each night. 

3. RESULTS 

3.1. Morphology 

To search for and characterize faint activity for 238P, we obtained multiple exposures of the 
comet each night. Nightly images were then shifted (using linear interpolation) to place the comet 
in the same location and then added together to create deep composite images (Fig. 1). These 
deep composite images chronicle a rise in activity for 238P over the course of these observations. 
While in principle, residual dust could remain from the comet's 2005 outburst, there is a clear 
transition from a nucleus-dominated morphology in July, August, and September (Fig. la-c) to a 
coma-dominated morphology in October and November (Fig. ld-e), the latter of which strongly 
resembles the morphology observed in 2005 (Hsieh et al. 2009b). As shown in Figure 2, these 
observations tightly constrain the time and orbital position of the onset of 238P's current active 
episode. This morphological transition is a strong indication that renewed dust production took 
place during this period, a key finding for ascertaining the origin of that dust production, and an 
issue we consider further in §3.2.1. 

3.2. Photometry 

3.2.1. Phase Function Analysis 

The variation of an object's brightness with solar phase angle (i.e., its phase function) is de- 
pendent on properties of that object's regolith, including albedo, particle size distributions, surface 
roughness, and porosity (c/. Helfenstein & Veverka 1989; Muinonen et al. 2002). Comparisons of 
phase function parameters found for different objects therefore provide a means for comparing the 
objects' surface properties themselves. 

Using data obtained in July and August 2010 when 238P was observed to be largely inactive, we 
obtain the best-fit phase function plotted in Figure 3a, defined by an absolute magnitude (at a = 0° 
and R = A = 1 AU) of H = 19.05±0.07 mag and slope parameter G = -0.03±0.10 in the IAU H, G 
system (corresponding to a linear phase-darkening coefficient of (3 ~ 0.045 ±0.007 mag deg -1 ). For 
this fit, we omit the inactive but extremely faint January 2007 observation from Hsieh et al. (2009b) 
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which we judged to be unreliable. For comparison, the slope parameters calculated for MBCs 133P 
and 176P are G = 0.04±0.05 and G = 0.26±0.05 (Hsieh et al. 2009a), respectively, and the range of 
slope parameters measured for other kilometer-scale Themis asteroids is —0.23 < G < 0.60 (Hsieh 
et al. 2008). This result therefore reinforces findings to date that aside from their cometary activity, 
MBCs exhibit no observable physical differences from dynamically similar but inert members of 
the background asteroid population. Using a geometric i?-band albedo of pr = 0.05, consistent 
with that of other MBCs (Hsieh et al. 2009a), we estimate an effective nucleus radius for 238P of 
r n ~ 0.4 km. 

Comparison of the remainder of our 2010 data to magnitudes predicted for an inactive 238P 
nucleus using the phase function derived here reveals a steadily increasing photometric excess for 
data from September through December 2010, corroborating the conclusion of our morphological 
analysis (§3.1) that 238P experienced renewed dust production, likely to be sublimation-driven (see 
§3.2.2), during this period. Plotting data (corrected to normalized photometry apertures) from 
both 2005-2007 and 2010 (Table 2) as a function of true anomaly (Fig. 3b), we find comparable 
photometric excesses in both data sets, with the excess observed in 2010 appearing poised to rise to 
the level observed in 2005. Further observations will be required when 238P reaches a comparable 
point in its orbit in late 2011 to confirm whether its current activity actually reaches a level 
comparable to its previous active episode. 



3.2.2. Sublimation Modeling 

To further analyze the photometric excesses detected in our data (§3.2.1), we employ a thermal 
sublimation model to estimate the grain flux, which can then be used to derive the estimated 
photometric enhancement expected from a sublimating surface (Meech et al. 1986). This model is 
capable of simulating multiple volatile components, but previous thermal modeling has shown that 
crystalline water ice is the only cometary volatile that is likely to be present on MBCs (Schorghofer 
2008; Prialnik & Rosenberg 2009). As such, for 238P, we only consider the action of H2O ice 
sublimation. 

For the purposes of the model, we assume that the comet's activity profile is identical during 
each apparition, permitting us to combine observations from 2005 and 2010 by parameterizing them 
by true anomaly. Due to the highly underconstrained nature of this modeling effort, we do not 
conduct a systematic search of the entire available parameter space, but instead seek to simply 
find a plausible sublimation scenario that reproduces our limited observations. Using standard 
assumptions for albedo and emissivity used for other comets and primitive asteroids, we find that 
our data fit plausibly well to a model in which activity turns on at v = 295° (corresponding to 
06 December 2004, or 21 July 2010), rising exponentially until reaching a maximum (scaled to 
match the maximum activity intensity) at v = 35° (24 November 2005, or 06 July 2011), and then 
immediately declining exponentially and ceasing at v = 120° (07 January 2007, or 08 August 2012). 
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Given the numerous underconstrained free parameters in our model, these activity initiation, 
peak, and cessation points certainly do not represent a unique solution and therefore do not nec- 
essarily reflect the precise reality of this comet's active behavior. They do, however, represent a 
plausible scenario in which 238P's activity is driven by the sublimation of water ice, accomplishing 
the primary goal of our modeling efforts. 

From our models, we find a likely active area equivalent to ~ 1% of 238P's total surface area, 
which for a r n ps 0.4 km body, represents an active area of ~ 2 x 10 4 m 2 , equivalent to a circular 
region with a radius of ~ 80 m. This radius is about an order of magnitude larger than the radius 
of the active area estimated for 133P (Hsieh et al. 2004). Along with 238P's smaller size and 
therefore lower escape velocity for emitted dust particles, the larger physical size of its active area 
is consistent with its dramatically more vigorous activity compared to 133P. 

Using the models, we also derive a maximum launchable dust grain size of ~ 10 mm, assuming 
a freely-sublimating surface. If sublimation is concentrated in vents, even larger grains could 
potentially be emitted. As such, long-lasting residual activity due to slow-dissipating large grains 
from 238P's previous active apparition in 2005 cannot be ruled out by the results of this modeling. 
Whether remnants of previous dust emission episodes are present or not, however, we already find 
strong evidence of renewed dust production in analyses discussed in §3.1 and §3.2.1. 

The delay in peak activity for 238P until ~ 35° past perihelion as implied by the best-fit 
model suggests either (1) the ice reservoir driving the observed cometary activity is located at 
some depth below the surface, delaying the onset of activity until a solar radiation-driven thermal 
wave of sufficient intensity is able to propagate to that depth, or (2) activity is modulated by 
seasonal effects as hypothesized for 133P (Hsieh et al. 2004) where maximum solar insolation of 
an isolated active site is primarily a function of an object's obliquity and solstice position, not 
heliocentric distance. If (1) is correct, we would expect all MBCs to exhibit peak activity shortly 
after perihelion. If (2) is correct, we would expect other MBCs to be equally likely to exhibit peak 
activity before perihelion as well. Thus far, all the known MBCs primarily exhibit post-perihelion 
activity, though better characterization of precise activation and cessation points, as well as the 
discovery and characterization of a larger sample MBCs, are required to determine whether (1) or 
(2) is the primary mechanism for MBC activity modulation. 

4. DISCUSSION 

Our finding of renewed dust production is strong confirmation that 238P is exhibiting genuine 
cometary activity driven by the sublimation of volatile ice. As concluded for 133P (Hsieh et al. 
2004), repeated activity is extremely unlikely to occur if 238P's observed dust emission is simply an 
ejecta cloud produced by an impactor striking inert asteroidal regolith. This scenario would require 
two separate collisions to occur on 238P within less than 5 years, when collisions causing similar 
behavior are not observed anywhere nearly as frequently on other nearby asteroids. The extraor- 
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dinarily coincidental timing required to produce two active episodes at nearly the same portions 
of 238P's orbit further decreases the plausibility of impact generation for its activity. A far more 
plausible explanation is that 238P's activity is driven by the sublimation of volatile ice, perhaps 
preserved in a subsurface reservoir still reachable (after a delay) by solar insolation, or contained 
within an isolated area exposed on the surface via collisional excavation and subsequently subjected 
to seasonal variations (either due to the object's obliquity, in direct analogy to Earth's seasons, or 
possibly also related to the local topography surrounding the active site) in solar radiation. 

Observations of 238P's activity in 2010 are thus far insufficient to evaluate the suggestion made 
by Hsieh et al. (2009b) that we had perhaps witnessed the immediate aftermath of a collisional 
activation of the MBC in 2005 and that subsequent outbursts may be much weaker. The obser- 
vations of 238P reported here were made at a different portion of its orbit than in 2005 and it is 
unclear whether the activity will increase or decline in intensity as the object progresses to that 
later portion of its orbit again. For this reason, further observations of 238P are highly encouraged 
when it again becomes observable from Earth in mid-2011. 

Even if activity levels during the current active episode remain comparable to those during 
2005, it is still possible that 238P was activated much more recently than less vigorously-active 
MBCs like 133P and 176P, except that its activity declines on a more gradual basis than from one 
orbit to the next. Long-term monitoring of 238P (and other vigorously-active MBCs such as P/2008 
Rl (Garradd) and P/2010 R2 (La Sagra)) will help clarify this behavior and likely give valuable 
insights into the scale of volatile reservoirs in active MBCs, and by extension, in compositionally 
similar but otherwise inactive main-belt asteroids, as well as their rate of depletion once excavated. 

We find that in comparison to 133P, whose periodic cometary activity appears to begin shortly 
before perihelion (Hsieh et al. 2010), 238P shows periodic activity beginning much earlier relative to 
perihelion. This discrepancy is consistent with the seasonal hypothesis proposed for the modulation 
of activity in 133P (and by extension, in other MBCs) (Hsieh et al. 2004), as it shows that MBC 
activity is not strictly correlated to heliocentric distance. This interpretation is not conclusive 
proof of this hypothesis, but does highlight the need for regular monitoring of other MBCs. Should 
MBCs consistently show activity only near perihelion, as most of the currently known MBCs do, 
this would be a sign that sublimation-driven cometary activity in these bodies is in fact strongly 
dependent on heliocentric distance, perhaps subject to moderate delays in the propagation of solar 
thermal waves introduced by insulating surface material. The collisional activation model for MBCs 
and correlated seasonal modulation hypothesis suggest, however, that other MBCs may exhibit 
activity well before perihelion, since under those hypotheses, activity modulation depends more on 
shadowing due to obliquity and topography than heliocentric distance (though some observational 
bias towards detections of active MBCs at smaller geocentric, and therefore heliocentric, distances 
is expected). 

HHH is supported by NASA through Hubble Fellowship grant HF-51274.01 awarded by the 
Space Telescope Science Institute, which is operated by the Association of Universities for Research 
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Table 1. Observation Log 



UT Date 


Tel. 


Moon a 


9 s h 


N c 


R d 


A c 


a f 








Active? j 


2010 Jul 07 


UH2.2 


N-5 


077 


16 


2.704 


1.821 


13.0 


291.8 


253.5 


251.0 


no 


2010 Jul 20 


UH2.2 


N+8 


(y/7 


23 


2.674 


1.709 


8.5 


294.8 


256.0 


251.4 


no 


2010 Aug 15 


SOAR 


N+5 


l'/O 


21 


2.616 


1.608 


2.6 


301.1 


55.0 


252.9 


no 


2010 Sep 03 


UH2.2 


N-5 


(y/7 


5 


2.576 


1.643 


10.7 


305.9 


69.9 


254.0 


no 


2010 Sep 04 


NTT 


N-4 


i'/i 


54 


2.574 


1.647 


11.0 


306.1 


70.1 


254.1 


yes 


2010 Sep 05 


NTT 


N-3 


078 


35 


2.572 


1.651 


11.4 


306.4 


70.3 


254.1 


yes 


2010 Oct 05 


Keck 


N-2 


(y/9 


8 


2.514 


1.869 


20.3 


314.3 


73.0 


254.4 


yes 


2010 Nov 25 


UH2.2 


N-11 


(y/8 


8 


2.433 


2.414 


23.5 


328.5 


71.3 


251.0 


yes 


2010 Dec 09 


UH2.2 


N+3 


0"8 


10 


2.416 


2.566 


22.5 


332.5 


70.3 


249.7 


yes 



a Phase of moon, in offset from new moon ("N") in days 
b Average seeing (FWHM) in arcsec 
c Number of images 
d Heliocentric distance (AU) 
e Geocentric distance (AU) 

f Solar phase angle (Sun-238P-Earth) in degrees 
g True anomaly in degrees 

h Position angle of the antisolar vector, as projected in the plane of the sky, in degrees east of north 

'Position angle of the negative velocity vector, as projected in the plane of the sky, in degrees east 
of north 

Ts visible activity detected? 
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Table 2. Photometric Analysis 



UT Date 




m Rtm i d b 


m R (l, l,a)° 


6m R d 




Mi 


2005 Nov 10 g 


19.28±0.05 





16.5±0.4 


2.7±0.4 


12±4 


10±4 


2005 Nov 19 g 


19.34±0.05 





16.6+.0.4 


2.9±0.4 


15±4 


13±4 


2005 Nov 20 s 


19.46±0.05 





16.7+.0.4 


2.8±0.4 


13±4 


11±4 


2005 Nov 218 


19.37±0.05 





16.6+.0.4 


3.0±0.4 


16±5 


14±4 


2005 Nov 22S 


19.28±0.05 





16.5+.0.4 


3.1+.0.4 


17±5 


15±4 


2005 Nov 26 s 


19.72±0.05 





16.9±0.4 


2.8±0.4 


13±4 


11±4 


2005 Dec 24§ 


20.12±0.03 




16.9+.0.4 


3.3±0.4 


21±6 


18±5 


2DD5 Dec 25§ 


on i fi+n D3 




i 7 n_|_n 4 


3.2+.0.4 


1 9+fi 


1 7+5 


2007 Jan 27 g 


24.9±0.4 




20.2±0.6 


-0.6±0.6 






2010 Jul 07 


23.61±0.10 


23.6±0.4 


20.1±0.4 


-0.1±0.4 






2010 Jul 20 


22.85±0.06 


22.9±0.4 


19.6±0.4 


0.2±0.4 






2010 Aug 15 


22.34±0.05 


22.5±0.1 


19.4±0.1 


0.0±0.1 






2010 Sep 03 


21.97±0.04 


22.0±0.4 


18.9±0.4 


1.0±0.4 


2.5±1.7 


2.2±1.5 


2010 Sep 04 


22.01±0.05 


22.3±0.2 


19.2±0.2 


0.7±0.2 


1.9±0.3 


1.6±0.3 


2010 Sep 05 


22.02±0.05 


22.3±0.2 


19.2±0.2 


0.7±0.2 


1.9±0.3 


1.6±0.3 


2010 Oct 05 


22.25±0.05 


22.3±0.4 


18.9±0.4 


1.4±0.4 


3.6±1.1 


3.1±1.0 


2010 Nov 25 


21.75±0.05 


21.8±0.4 


17.9±0.4 


2.5±0.4 


10±3 


8.7±2.6 


2010 Dec 09 


21.86±0.07 


21.9±0.4 


17.9±0.4 


2.5±0.4 


10±3 


8.7±2.6 



a Mean .R-band magnitude (averaged in flux space), measured in 4'/0 radius apertures 

b Estimated midpoint of i?-band lightcurve, assuming peak-to-trough photometric 
range of Am ~ 0.8 mag 

c Reduced midpoint R-band magnitude (normalized to R = A = 1 AU) 

d Deviation between measured m R (l,l,a) and predicted magnitude from best-fit 
phase function based on data from July- August 2010 when 238P was observed to be 
inactive 

c Inferred ratio of scattering surface area of dust to nucleus scattering surface area 

f Estimated dust mass, in 10 4 kg, assuming 10 //m-radius grains and p=1300 kg m~ 3 , 
as for 133P (Hsieh et al. 2010), contained within 4'/0 photometry apertures 

sProm Hsieh et al. (2009b). 



Fig. 1. — Composite i?-band images of 238P (at the center of each panel) constructed from data obtained 
on (a) UT 2010 July 7 (8100 s of effective exposure time on the UH 2.2 m), (b) UT 2010 July 20 (19530 s on 
the UH 2.2 m), (c) UT 2010 August 14 (4050 s on SOAR), (d) UT 2010 September 03 and 04 (23100 s on 
the NTT), (e) UT 2010 October 05 (840 s on Keck), (f) UT 2010 November 25 (4800 s on the UH 2.2 m), 
and (g) UT 2010 December 09 (4500 s on the UH 2.2 m). All panels are 6070 x 60'.'0 in size, with North 
(N), East (E), the antisolar direction (—0), and the negative heliocentric velocity vector (— v), as projected 
on the sky, marked. 
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Fig. 2. — Orbital position plot for 238P observations detailed in Table 1 and in Hsieh et al. (2009b). The 
Sun is shown at the center as a solid dot, with the orbits of Mercury, Venus, Earth, Mars, 238P and Jupiter 
(from the center of the plot outwards) shown as black lines. Green squares mark positions where 238P was 
observed to be active in 2005 and red circles mark active positions in 2010. Open squares mark positions 
where no activity was detected between 2007 and 2010. Perihelion (P) and aphelion (A) positions are marked 
with crosses. Observations plotted are from (a) 2005 October 24 (Read et al. 2005), (b) 2005 November 
10, (c) 2005 November 19-22, (d) 2005 December 24-25, (e) 2007 January 27, (f) 2010 July 07-20, (g) 2010 
August 15, (h) 2010 September 03-05, (i) 2010 October 05, (j) 2010 November 25, and (k) 2010 December 
09, where (b)-(e) are from Hsieh et al. (2009b) and (f)-(k) are from this work. 



22 



5 10 15 20 

Solar phase angle (deg) 




-100 -50 50 

True anomaly (deg) 



100 



Fig. 3. — (a) Best-fit phase function (solid line) for reduced i?-band magnitude data (normalized to R = 
A = 1 AU) from July-August 2010 (open squares) when no activity was detected for 238P (January 2007 
photometry is also plotted for reference) with the range of uncertainty due to phase function parameter 
uncertainties shown by dotted lines. Data from November-December 2005 (green squares; Hsieh et al. 
2009b) and September-December 2010 (red circles) when 238P was active are overplotted. (b) Plot of results 
of sublimation modeling where photometry points from (a) are normalized to a — 0° and to fixed photometry 
apertures of 5'.'0 at R = A = 1 AU (assuming a r _1 surface brightness profile during active periods) and 
replotted, the baseline magnitude of the nucleus is marked by a thick black line (where the phase function is 
approximated as linear), and the brightness inside a 5'.'0 photometry aperture predicted by our sublimation 
model is marked by a thick blue line. Error bars in both panels indicate the range of uncertainty due to 
unknown rotational phases and standard photometric uncertainty, assuming a peak-to-trough photometric 
range of Am ~ 0.8 mag. Given the large uncertainty associated with our January 2007 photometry (at 
v = 123°), it is not explicitly used to fit our model, though we do assume the comet to be inactive at that 
point. 



